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Abstract. The kinetic theory is used to study the evolution of the self-gravitating 
£SJ . disk of planetesimals. The effects of frequent collisions between planetesimals are 

taken into account by using a Krook integral in the Boltzmann kinetic equation. 
It is shown that as a result of an aperiodic collision-dissipative instability of small 
gravity disturbances the disk is subdivided into numerous dense fragments. These 
can eventually condense into the planetary sequence. 
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Solar system formation is thought to start with dust particles set- 
tling to the central plane of a nebula to form a thin dust layer around 
the equatorial plane. During the early evolution of such a rapidly rotat- 
ing disk it is believed that the dust particles coagulate into numerous 
kilometer-sized rocky bodies (planetesimals). See Taylor (1992) as a 
review of the problem. 

In a swarm of planetesimals direct physical collisions inevitably be- 
come an important factor (Taylor 1992). One can suggest that planets 
(Mercury, Venus, . . ., Neptune) accreted subsequently from a hierarchy 
of colliding planetesimals. (The combination of low mass with a highly 
inclined and eccentric orbit is a major reason for not according Pluto 
planetary status: observations point to the orbit of Neptune as the true 
outer boundary of the planetary system.) 

Our principal idea is to regard the formation of the planetary system 
as a possible last stage in the formation of the self-gravitating, highly 
flattened solar nebula with frequent, almost elastic collisions between 
rocky planetesimals. We argue that a collision-dissipative collective in- 
stability of small-amplitude gravity perturbations developing in such a 
disk leads to the formation of an arrangement of dense aggregates. We 
speculate that subsequent substantial gravitational interaction between 
these dense aggregates circling the primordial sun would result in the 
formation of the planetary sequence. When cores of the giant planets 
reach a critical mass (~ 10 masses of the Earth) they begin to accrete 
the interplanetary gas. 

The collision motion of an ensemble of identical planetesimals in the 
plane, in the frame of reference rotating with angular velocity $7, can 
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be described by the Boltzmann kinetic equation 

dt df df dv \dt J coll ' 

where f(f,v,t) is the phase-space distribution function of planetesi- 
mals, $(r, t) is the total gravitational potential, and (df/dt) co \\ is the 
so-called collision integral which takes into account effects due to the 
discrete-point nature of the gravitational charges and defines the change 
of / arising from ordinary interparticle collisions. There are only elastic 
physical collisions, and momentum is conserved in collisions. There is no 
correlation in motion between the colliding species, that is, Boltzmann's 
hypothesis of molecular chaos is adopted. 

The Poisson equation in a suitable form is V 2 <£ = A / wGa5(z), where 
a(f, t) is the surface mass density. The Boltzmann and Poisson equa- 
tions, with the distribution function of planetesimals and appropriate 
boundary conditions, give a complete description of the problem for 
infinitesimally thin disk modes in a self-consistent field problem. In 
the present study, this system of equations is simplified by consider- 
ing the lowest WKB approximation (in the limit v c 3> f2, where v c 
is the frequency of collisions between planetesimals); this is accurate 
for short wave perturbations only, but qualitatively correct even for 
perturbations with a longer wavelength, of the order of the system 
radius. 

As is known, the Boltzmann equation is nearly intractable because 
of the complicated collision integral. This integral can be approximated 
in various ways. In this work we use the simplified kinetic model when 
the exact integral (df /dt) co \\ is replaced by a phenomenological term 
in the form of a Krook model (Lifshitz Sz Pitaevskii 1981; Griv Sz Peter 
1996; see also Griv et al. 1999). 

Perturbations in the gravitational field cause perturbations to the 
planetesimal distribution function. In the linear approximation, one can 
therefore write f(f,v,t) = fo(r,v) + fi(f,v,t), where |/i| < f and fi 
is a function rapidly oscillating in space and time. We consider zeroth 
order disk inhomogeneity in the r-direction only. Initially the disk is 
in equilibrium, dfo/dt = 0. The function /q describes the rotating 
"background" against which small perturbations develop. If an initial 
perturbation grows, the system is called unstable. 

As a result of this study, the dispersion relation connecting the 
Doppler-shifted wavefrequency cj* and the wavenumber k of pertur- 
bations is found: 
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where <Jo(r) is the equilibrium surface density, c is the dispersion of 
random velocities of planetesimals ("temperature"), g = 4(7 + v c ), and 
erf(/3) denotes the error function. The wavefrequency is represented in 
the form u = ij, where 7 is real and positive (cf. Binney & Tremaine 
1987, p. 292). 

A numerical solution of the dispersion relation is shown in Figure 1. 
As follows from this solution, the dispersion relation gives the condition 
for an aperiodic collision-dissipative instability (7 > 0): 

A > A cr it ~ c 2 /Go"o and A = 2ir/k, 

which is just the familiar Jeans criterion for a nonrotating medium and 
A cr it is the ordinary Jeans length (Binney & Tremaine 1987). Frequent 
collisions thus remove the rotational stabilization in a flat system. Thus, 
the theory predicts that a gravitating disk with frequent interparticle 
collisions is aperiodically unstable for perturbations with wavelengths 
greater than the Jeans length. 

At the wavelength Gao / c 2 the growth rate of the instability is 
maximum (Fig. 1). It means that of all harmonics of initial perturba- 
tion, one perturbation with the maximum of the growth rate and with 
k* will be formed asymptotically in time. The numerous condensations 
with typical dimensions and distances between them A* = that 
arise will remain localized in space and grow, because the instability 
will be an aperiodic. 
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Figure 1. The dispersion relation for reasonable parameters of a two-dimensional 
system of frequently colliding planetesimals. 
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